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ABSTRACT 



Aims. We re-examine the possibility that astrophysical jet collimation may arise from the geometry of rotating black holes and the 
presence of high-energy particles resulting from a Penrose process, without the help of magnetic fields. 

Methods. Our analysis uses the Weyl coordinates, which are revealed better adapted to the desired shape of the jets. We numerically 
integrate the 2D-geodesics equations. 

Results. We give a detailed study of these geodesies and give several numerical examples. Among them are a set of perfectly colli- 
mated geodesies with asymptotes p = pi parallel to the z-axis, with pi only depending on the ratios and ^ , where a and M are 
the parameters of the Kerr black hole, E the particle energy and Q the Carter's constant. 

Key words, black hole physics; relativity; galaxies: jets; accretion; accretion discs 



1. Introduction 

It has long been speculated that a single mechanism might be 
at work in the production and collimation of various very ener- 
getic observed jets, such as those in gamma ray bursts (GRB) 
KPiranet aL, 0T)1 |(Sheth et al.,"03)l |(Fargion, 03)1 and jets 



ejected from active galactic nuclei (AGN) (Sauty et al., 02) 
and from microquasars (Mirabel & Rodriguez, 94) 



[(Mirabel & Rodriguez, 9 9) Here we limit ourselves to jets pro 
duced by a black hole (BH) type core. The most often invoked 
process is the Blandford-Znajek (Blandford & Znajek, 77) or 



some closely similar mechanism (e.g. (Punsly & Coroniti, 90a) 
(Punsly & Coroniti, 90b) (Punsly, 01)|) in the framework of 



magnetohydrodynamics, always requiring a magnetic field. 
However such mechanisms are limited to charged particles, and 
would be inefficient for neutral particles (neutrons, neutrinos 
and photons), which are cuiTently the presumed antecedents 
of very thin and long duration GRB (Fargion, 03) Moreover, 
even for charged particles, some questions persist (see for 
instance the conclusion of |(Williams, 04)[ ). Finally, while the 
observations of synchrotron radiation prove the presence of 
magnetic fields, they do not prove that those fields alone cause 
the collimation: magnetic mechanisms may be only a part of a 
more unified mechanism for explaining the origin and collima- 
tion of powerful jets, see ,(Livio, 99)| p. 234 and section 5, and, 
in particular, for collimation of jets from AGN to subparsec 
scales, see |(de Felice & Zanotti, 003| 

Considering this background, it is worthwhile looking for 
other types of model to explain the origin and structure of 
jets. Other models based on a purely general relativistic ori- 
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gin for jets have been considered. A simple model was ob- 
tained by |(Opher et al., 9 6)'by assuming the centres of galaxies 
are described by a cylindrical rotating dust. That paper showed 
that confinement occurs in the radial motion of test particles 
while the particles are accelerated in the axial direction thus 
producing jets. Another relativistic model was put forward in 
|(Herrera & Santos, 07)| This showed that the sign of the proper 
acceleration of test particles near the axis of symmetry of quasi- 
spherical objects and close to the horizon can change. Such an 
outward acceleration, that can be very big, might cause the pro- 
duction of jets. 

However, these models show a powerful gravitational effect 
of repulsion only near the axis, and are built in the framework of 
axisymmetric stationary metrics which do not have an asymp- 
totic behaviour compatible with possible far away observations. 
So we want to explore the more realistic rotating black hole, i.e. 
Kerr, metrics instead. 

We thus address here the issue of whether it is pos- 
sible, at least in principle (i.e. theoretically) to obtain a 
very energetic and perfectly collimated jet in a Kerr black 
hole spacetime without making use of magnetic fields. Other 
aut hors (see |(Bicak et al. 93)] [(Williams, 95)] [(Williams, 04")] 
(de Felice & Carlotto, 97) and references therein) have made 
related studies to which we refer below. Most such authors 
agree that the strong gravitational field generated by rotating 
BHs is essential to understanding the origin of jets, or more 
precisely that the jet originates from a Penrose-like process 
(Pe nrose, 69)[[(WiUiams, 04)| in the ergosphere of the BH; colli- 
mation may also arise from the gravitational field and that is the 
main topic in this paper. 

Our work can therefore be considered as covering the whole 
class of models in which particles coming from the ergosphere 
form a jet collimated by the geometry. Although a complete 
model of an individual jet would require use of detailed models 
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of particle interactions inside the ergosphere, such as that given 
by (Will iams, 04)| we show that thin and very long and energetic 
jets, with some generic features, can be produced in this way. In 
particular the presence of a characteristic radius, of the size of 
the ergosphere, around which one would find the most energetic 
particles, might be observationally testable. 

From a strictly general relativistic point of view, test particles 
in vacuum (here, a Kerr spacetime) follow geodesies; this applies 
to both charged and uncharged particles, although, of course, 
in an electrovacuum spacetime, such as Kerr-Newman, charged 
particles would follow accelerated trajectories, not geodesies. 
Thus, in Kerr fields, what produces an eventual collimation for 
test particles, or not, is the form of the resulting geodesies. 
Hence we discuss here the possibilities of forming an outgoing 
jet of collimated geodesies followed by particles arising from 
a Penrose-like process inside the ergosphere of a Kerr BH. We 
show that it is possible in principle to obtain such a jet from 
a purely gravitational model, but it would require the "Penrose 
process" to produce a suitable, and rather special, distribution of 
outgoing particles. 

The model is based on the following considerations. 

Most studies of geodesies, e.g. |(Chandrasekhar, 1983)| em- 
ploy generalized spherical, i.e. Boyer-Lindquist, coordinates. 
We transform to Weyl coordinates, which are generalized cylin- 
drical coordinates, and are more appropriate, as we shall see, for 
interpreting the collimated jets. 

We consider test particles moving in the axisymmetric sta- 
tionary gravitational field produced by the Kerr spacetime, 
whose geodesic equations, as projected into a meridional plane, 
are known |(Chandrasekhar, 1983)| Our study is restricted to 
massive test particles, moving on timelike geodesies, but of 
course massless test particles on null geodesies could be the 
subject of a similar study. (Incidentally the compendium of 



(Sharp, 79) shows that analytic studies of general timelike 



geodesies have been much less frequent than detailed studies of 
more restricted problems.) 

For particles outgoing from the ergosphere of the Kerr BH 
we examine their asymptotic behaviour. Among the geodesic 
particles incoming to the ergosphere, we discuss only the ones 
coming from infinity parallel to the equatorial plane, because 
these are in practice the particles stemming from the accretion 
disk. We show that only those with a small impact parameter are 
of high enough energy to provide energetic outgoing particles. 

In the ergosphere, a Penrose-like process can occur In the 
original Penrose process, an incoming particle decays into two 
parts inside the ergosphere. It could also decay into more than 
two parts, or undergo a collision with another particle in this 
region, or give rise to pair creation (e",e^) from incident pho- 
tons which would follow null geodesies. The different possi- 
ble cases do not affect our considerations, and that is why 
we do not study them here, although the distribution func- 
tion of outgoing particles would be required in a more de- 
tailed model of the type discussed, in particular to explain why 
only particles with low angular momentum and not diverging 
from the rotation axis are produced. For detailed studies see 
[(Wilhams, 95)i |(Williams, 4), (Pira n & Shaham, 77)| After a 
decay, one (or more) of the particles produced crosses the event 
horizon and irreversibly plunges into the BH, while a second 
particle arising from the decay can be ejected out of the ergo- 
sphere following a geodesic towards infinity. This outgoing par- 
ticle could be ejected so that asymptotically it runs parallel to the 
axis of symmetry, but we do not discuss only such particles. 

In our model there is no appeal to electromagnetic forces to 
explain the ejection or the collimation of jets, though the parti- 



cles therein may themselves be charged. The gravitational field 
suffices, in the case of strong fields in general relativity, which is 
the case near the Kerr BH, provided the ergosphere produces par- 
ticles of appropriate energy and initial velocity. The gravitomag- 
netic part of the gravitational field then provides the collimation. 
Hence, our model is, in this respect, simpler than the standard 
model of (Blandford & Znajek, 77) and is in accordance with 
the analysis given in .(Williams, 04)| 

The paper starts with a study of Kerr geodesies in Weyl co- 
ordinates in section |2j the next section studies the asymptotic 
behaviour of geodesies of outgoing particles with L, - 0; sec- 
tion|4]analyses incoming particles stemming from the accretion; 
a sample Penrose process and the plotting of geodesies are pre- 
sented in section |5j and finally we discuss in section |6] the sig- 
nificance of our results for jets. In the conclusion, we succinctly 
summarize our main results and evoke some perspectives. 

2. Kerr geodesies 

We start from the projection in a meridional plane (p - constant 
of the Kerr geodesies in Boyer-Lindquist spherical coordinates 
r, 9 and S. The metric is 



ds^ = (r^fl^cos^e) 



sin^e 



- 2Mr + a^- 



(r^ + d?- cos^ 8) 
{f - 2Mf + a^) 
(f^ + d?- cos^ 6) 



[aAt - (r 
(dt - a : 



(1) 



where M and Ma are, respectively, the mass and the angular mo- 
mentum of the source, and we have taken units such that c - I = 
G where G is Newton's constant of gravitation. The 'radial' co- 
ordinate in ([T]i has been named ? because it is more convenient 
for us to use the rescaled coordinate r = r/M , (O'Neill, 95)| The 
projected timelike geodesic equations are then 



(fl^r -I- a%r + arr + a\r + an) 
Mr - , 



M^e^ 



r2 + (§) cos2 
b4 cos** 6 + b2 cos^ + bo 



(1 -cos^fl) 



r2 + (§)'cos2 



(2) 



(3) 



with coeflicients 
a^Q 



flo = - 



M4 ' 



fl2 - 

fl3 - 
fl4 - 

and 

bo - 
bi - 
b4 - 



'[a\E^--l)-Ll-Q]., 



M2 

2, 



Q_ 



1)-Ll-Q]^a2, 



1); 



(4) 

(5) 

(6) 

(7) 
(8) 

(9) 
(10) 
(11) 
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where the dot stands for differentiation with respect to an afiine 
parameter and E, L, and Q are constants. Here Chandrasekhar's 
6\ has been set to 1, its value for timelike curves. Assuming 
that the affine parameter is proper time t along the geodesies, 
then these equations implicitly assume a unit mass for the 
test particle Q so that E and L- have the usual significance 
of total energy and angular momentum about the z-axis, and 
Q is the corresponding Carter constant (which, as described 
in [(Hughson et al., 72)| for example, arises from a Killing ten- 
sor of the metric, while E and L, arise from Killing vectors). 
With this understanding, E, L,, and Q have the dimensions 
of Mass, Mass^ and Mass^ respectively, in geometrized units, 
while 6i, though 1 numerically, has dimensions Mass^, as do 
all the a, and bi. In this paper we consider only particles on 
unbound geodesies with E > I. (For the conditions for exis- 
tence of a turning point, giving bound geodesies, which are re- 
lated to the parameter values for associated circular orbits, see 
KChandrasekhar, 1983)| [(Williams, 95)||(Wimams, 04)| ) 

The dimensionless Weyl cylindrical coordinates, in multiples 
of geometrical units of mass M, are given by 



P 



[{r-lf -Af'^sinO, z = (r - l)cos6l, 



where 



From (fT2l l we have the inverse transformation 



r = a + 1 
sin 6 = 



— COS0=— , 

(a2-A)i/2 a 



with 



(12) 

(13) 

(14) 
(15) 

(16) 



Here we have assumed A > 0, and taken the root of the second 
degree equation obtained from ( fTSl l for the function a^(p, z) that 
allows the extreme black hole limit A = 0. The other root, in this 
limit, is the constant a = 0. 

The equation ( fTSI l shows that in the (p, z) plane the curves of 
constant a (constant r) are ellipses with semi-major axis a and 
eccentricity e - y/A/a: for large a, these approximate circles. 
Note that p - consists of the rotation axis 6 - or n together 
with the ergosphere surface. 

Now, with ( fT4l i and ( fTSl ) we can write the geodesies (|2| and 
(O in terms of p and z coordinates, producing the following au- 
tonomous system of first order equations 



Pa^p _^ S(a^-A)z 



Mp = 



ap 



(„.i,V.(±) 



Mz = iPz-S)a 



(a+1) 



where 

P - e\ [04(0' + 1)"^ + a^{a + if" 

2 l'/^ 

-l-fl2(a-l-l) -H fli(Q' + 1) -H flo : 



S = 6162(^42'* + ^aaV + boa'*) 



.4^1/2 



(17) 



(18) 



(19) 



(20) 



' An alternative interpretation is to assume that for a particle of mass 
m, the affine parameter r/m has been used |(WiIkins, 72)||(Williams, 95)| 



and e, - ±1 for i - 1,2: ei indicates whether the geodesic is 
incoming or outgoing in r (i.e. the sign of r), while ei £2 indicates 
whether 9 is increasing or decreasing. Note we always mean the 
non-negative square roots to be taken. 

The ratio between the first order differential equations ( fTTI ) 
and (fTSl l yields the special characteristic equation of this system 
of equations 



dz ^ {\P\z-e2\S\)(a^-A)a^p 
dp \P\a'^p^ + £2\S\(a^-A)h' 



(21) 



We restrict our study to the quadrant p > and z > in 
the projected meridional plane (orbits in fact spiral round in (p 
in general: this information is contained in the conserved L^). 
The results for the other three quadrants will follow by sym- 
metry, although this symmetry does not imply that individual 
geodesies are symmetric with respect to the equatorial plane. 
From numerical solutions of the geodesies we obtained asym- 
metrical geodesies, confirming the analysis in [(Williams, 95)| 
[(Williams, 04)| Geodesies can also cross the polar axis, which 
would be represented by a reflection from p = back into the 
quadrant. 

Geodesies going to or coming from the expected accretion 
disk would, if the disk were thin, go to or from values of p much 
larger than z- In this limit (p^ z and p » VA) , we have 



a = p(l + 0(p-2)), 
\P\ = V^p2(l + k/p + 0{p-^)), 
2E^ - 1 

where k ■ 



£2 - 1 ' 

15 1 = ^/h)p\l + 0(fi-^)), 



and thus 



where 



dz 
dp 



Z - £2Zl + Zk/p 

P(l+k/p) 



+ 0{p-\ 



(bo 
Zl = — 
\a4 



1/2 



1 

M 



Q 



£2-1 



1/2 



(22) 
(23) 

(24) 

(25) 

(26) 
(27) 



and we have to assume <3 > to obtain a real zi (the form of (O 
makes it obvious that geodesies with <3 < are bounded away 
from the equatorial plane cos = 0, though those with very small 
\Q\ could still satisfy p » z). 

The truncated series development of dz/dp now yields 



dz 1 / 62^1 
T ~ - \Z- £2Zl + 

dp p\ p 



+ 0{p-\ 



(28) 



If 62 = - 1 the curve crosses the equatorial plane and gives simi- 
lar asymptotic behaviour in the second quadrant, so we can take 

62 = 1. 

A thicker accretion disk would absorb or release particles on 
geodesies with larger values of zip, which might include parti- 
cles with (3 < 0. 

Geodesies in an axial jet would have z ^ p- For this limit, 
we first observe that from ( [20l i we have 

S ^£2 [(bo + b2 + b4)z' + (2bo + b2)p'z'-f^ + 0(z-'), (29) 
where 

bo + b2+ b4 



2bo + b2 = -^[a\E^-l)-Ll+Q]. 



(30) 
(31) 



4 



J. Gariel et al.: Kerr Geodesies, the Penrose Process and Jet Collimation by a Black Hole 



Hence in this limit S is well defined and real for indefi- 
nitely small p/z only for L. - 0. The geodesies obeying this 
restriction, imposed after similar reasoning, were studied by 
|(Bicak et al., 93)| but in BL or Kerr-Schild coordinates. Here we 
re-examine these geodesies in the more revealing cylindrical co- 
ordinates. 

Before doing so, we may note that in contrast to geodesies 
with L- + 0, geodesies with > I and Lr -0 may lie arbitrarily 
close to the polar axis |(Carter, 1968)] For L, + 0, the value of 5 ^ 
at the axis is -Z'l2. < which is not allowed and thus there is 
some upper bound 0o on Q. The value of 5~ at = 7r/2 is hco?' so 
if £2 < there is also a lower bound Q\ on Q. 

3. Geodesies witli - 

We shall discuss unbounded {E^ > 1) outgoing geodesies. 
Corresponding incoming geodesies will follow the same curves 
in the opposite direction. 

For L- -Q, S- factorizes as 



5^ = (a" - z0«3a + fl (£ - l)z')/M- 



Hence S can only be zero at the symmetry axis, where cos 6-1, 
a = z andp = 0, or, if (3 < 0, at some z/a = (\Q\/a^[E^-l]y'^ = 
cos 0[ , say. Correspondingly 6 = only at the axis, at = 0i if 
Q <Q, and at a — > oo. 

Thus for Lj = and (3 < 0, geodesies which initially have 
6 < will become asymptotic to = 6i. The angle may be 
narrow if 

a\E^-l)-\Q\«:\Q\, (33) 

and then 0i « 1. Such geodesies may provide a conical jet, as 
discussed later. 

Our other polynomial, P^, can be written as 



4 ^ a 2 2a 



P^ = (£^ + — rr^ + — rr 



a 



Afi Afi 
2 Q 



(34) 



H-2r' + 2—r -lr-2r+ — 



M2 



^From this form it easily follows that any unbound geodesic 
(£^ > 1) with Lr = has at most one turning point in r (i.e. value 
such that /• = 0) and this, if it exists, lies inside the horizon (and 
a fortiori inside the ergosphere) (Stewart & Walker, 74) The ar- 
gument is very simple. If E^ > \ and (3 < 0, then P'^ is strictly 
positive for all r > 0. If (3 > 0, is negative at r = but positive 
at the outer black hole horizon (where -2r+a^/M^ = 0), so its 
one zero lies inside the black hole. This implies that unbounded 
outgoing geodesies followed by particles with L, = must come 
from the ergosphere. Correspondingly, geodesies incoming from 
infinity with L, = will fall into the ergosphere. 

Although there are no turning points of r, one can have turn- 
ing points of p, if £2 = -1- Such turning points are solutions of 
the equation 



Dip, z) = \P\ay ■ 



■A)\ 







(35) 



where D(p, z) is the denominator of (ISTT i with £2 = -1. At 
each of these turning points (p2, Z2), dz/dp — > 00, which means 
that the geodesies have a vertical tangent (parallel to the z- 
axis). Before reaching the turning point, these geodesies have 
dp/dz > 0, and, at any z, dz/dp > dz/dplo, where dz/dp|o is 
the slope of the curve (|35] |. and afterwards they have dp/dz < 0, 
implying that they subsequently cross the axis. They will then 
cross the curve ( l35T l again but from above in the (p, z) plane and 



hence with dz/dp < dz/dp|o, and afterwards stay in the region 
outside ( |35] |. 

For outgoing geodesies outside (l35T l which reach points at 
large z and p (» VA), then unless the ratio of z to p is very large 
(the case which we discuss next) or very small, approximating 
( l2n i gives dz/dp « z/p, so all such geodesies approximate p = 
Cz for suitable C, regardless of the sign of €2- 

In the limit z^ p and z » VA, 



a = zil+Oiz'^)), 
\P\ = y/^z^(l+k/z+0(z-^)), 

\S\ = ^J2bo + b2 pz(l + 0(z-^)), 

so the equation (1211) can be approximated by 

dz z(l+k/z) ^^/l 
dp p{\ +klz) + e2P\ \pz 



where 



(32) Pi = 



2/7(1 + b2 

fl4 



1/2 



= Pf 



1 + 



Q 



fl2(£2_ 1) 



1/2 



(36) 
(37) 

(38) 



(39) 



(40) 



and pe = a/M. Here pi is real if a^(E^ - I) + Q > but we see 
from ( l32b that for S to be real near the axis, this condition must 
be satisfied. 

In Figure [1] we show a plot of the values of dz/dp, using 
(l211 i. for £2 = -1, with the parameters M = 1, a = 1/2, £ = lO'*, 
<3 = -9 X 10^. The only asymptotes are parallel to the z axis at 
p = pi as expected from (l40l i. 




Fig. 1. Plot of the surface dz/dp = /(p, z) given by equation ( |2TI ) 
for an outgoing particle with the parameters £ = 10'*, L; = 
and (3 = -9 X 10^ for a black hole with parameters M = 1 and 
a = 0.5, in the case where £2 = -1. The points where dz/dp — > 
±00 correspond to the asymptotes given by the equation (l40l i. 
p=pi =0.49991. 

We also plot in Figure |2] a set of such outgoing geodesies 
obeying ( |2TI ). for the same values of the parameters of the BH 
(fl = 1/2,M = 1) and of the particle (L- = 0,(3 = -2.2xl0^£' = 
2.10^ so pi = 0.441588), but with diff'erent initial values of the 
position. The set of turning points of these geodesies is the curve 
defined by equation ( |35] |. For the rightmost of these geodesies, 
the numerical integration was also continued back towards the 
ergosphere as far as p = lO""*, z = 0.843407 < Va = 0.866025. 



J. Gariel et al.: Kerr Geodesies, the Penrose Process and Jet Collimation by a Black Hole 



5 




Fig. 2. Plots of geodesies obeying the equation (I2TI 1. showing the 
turning points. From left to right these curves start at po = 0.07 
and zo = 0.98, 0.95, 0.93, 0.92, 0.91, 0.9 and 0.89935501 

To confirm the picture obtained from these numerical exper- 
iments, one can show, without assuming z » p, the existence of 
exactly one zero of D on any curve r - constant, < 6 < n/l, 
so that the conclusion that a geodesic has at most one turning 
point in p is not an artefact of the approximation at large z- The 
argument is as follows. 

Along an r = constant curve, |P| and a are constant, p - 
posin0 and z = o-cos^, where po is a constant (related to r). 
Then 



D = IPla'^plsin^ - \S\(a^ -Af a cos 
where, defining 



F = (bo + \b4\cos^0y'^ 

= [(Q + a^(E^ - l)cos^ 0)]^''-/M, 



(41) 



(42) 



from ( |32] ) we have \S \ - a^F sin 0. For real S we need F > Q 
and as decreases, F increases. 

At 6* = 7r/2, D > 0, while at small 0,D<O. Hence there is at 
least one zero of D. Let the largest one be at 6* = 0o say. On the 
r - constant curve, we will then have 



D = a\a^ - AY 



sin0 
sin 0Q 



{F()Sm0cos0o - F cos0sm0o). (43) 



Here we have used D = at = 6*0 to substitute for the constant 
|P|aVo ^ii^^ ^0 in terms of a and 0o. As decreases from 0o, sin 
decreases, cos increases and F increases. Thus the combination 
(F() sin cos 00 - F cos sin ^o) becomes and stays negative for 
all < 00. Thus D < for all < 0o, though D approaches 0, 
due to the further factor sin 0, as — > 0. This implies that the 
only points on r = constant such that D = are at = and 
0^00- 

For large z we see from (39[ that the turning points lie ap- 
proximately on a curve p - piz/(z + k) or z = -kp/(p - pi). 
Actually, the differential equation for large z, if we drop the l/z^ 
terms, has an analytic solution 

Ckz ^kp + pizln[z/(z + k)], (44) 

where C is a constant of integration, so 

p = piz/feln(l + k/z) + Cz^ px+Cz... (45) 

as z —> 00. Keeping the next order terms in 1/z would be incon- 
sistent with the terms dropped during the derivation. (Similarly, 



there is an analytic solution for the approximate equation (l26b 
at large p which gives z - Zipln(l -1- p/k)/k + cp with similar 
interpretation.) 

^From ( |45] ), either (a) p/z is approximately constant or (b) 
p — » pi. In case (a), we note that for consistency of the ap- 
proximation z » p, C must be small, although the conclusion is 
the same as was reached above merely with the assumption that 
both z and p are » a/A. In case (b), we have a limit-outgoing 
geodesic for which p < pi at all points and as z — > co, p ^ pi. 
This limit is obtained since the turning point for p has Z2 — > c» 
when p2 — > pi. We can see from (l44l i that the coordinate Z2 of 
this turning point tends to infinity like Z3 = -kpi /(p2 - Pi )• The 
geodesies asymptotic to pi would provide a perfectly coUimated 
jet parallel to z. 

One might think (and we initially thought) that there also 
existed geodesies eventually tending to the same asymptote but 
approaching it from the right in the (p, z) plane (for example, 
directly from the accretion disk, or coming from the ergosphere 
but with a turning point p2 > pi). However, such geodesies do 
not exist, since they require that dz/dp < in the limit z » p 
and for p > pi , contradicting ( [39] l which implies dz/dp > 0. This 
is entirely in agreement with the results of Stewart and Walker 
[(Stewart & Walker, 74)[ 

The geodesies in p > pi may asymptote to any ratio z/p, 
from ( |39] |. Moreover, geodesies which do turn in p then cross 
the axis, cannot cross the curve D = from below again, and 
so cross it from above and also asymptotically have some fixed 
ratio z/p. 

For astrophysical applications, it may be important to write 
the results in the normal units of length and time. We have, from 
([TtT i. that asymptotically for outgoing particles in z » p, p — » 0, 
f — > £ , z > (for -62 = 1 = ei) and (fTSl l is given by 



Mz ~ V«4; 



(46) 



hence, restoring normal units of length and time and taking a par- 
ticle of mass m, the asymptotic value v of the speed of outgoing 
particles is given by 



mMz '■ 



1/2 



my- 



where we have used (O and 



r = 



-1/2 



(47) 



(48) 



is the Lorentz factor. Hence, asymptotically, the speed of the par- 
ticle is 

■• (49) 



V - c- 



= c 1- 



2 4\l/2 

m c * 



£2 



which is ultrarelativistic if £ » y/mc-. From (l47l l we have that 
asymptotically limit-outgoing particles have an uniform motion 
parallel to the z axis. 

4. Incoming particles 

We describe as "incoming particles" the particles, with param- 
eters E', L'., and Q' , coming into the ergosphere following un- 
bound geodesies and having a turning point in z (i.e. such that z 
- 0). Such turning points {p4, za] are defined as solutions of the 
equation 

A?2(P4,Z4) = (50) 
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where 



N2^Pz-S 



(51) 



is the relevant factor in the numerator of the right side of ( 1211 1. 
As remarked earlier we need only consider = 1 . For large p 
these turning points are approximately at z = zi(l - k/p), from 
(|28] ). The set of these turning points {p4, 14} forms a curve of the 
points of each geodesic with horizontal tangent (i.e. parallel to 
the p-axis). Each geodesic has a positive slope dz/dp > for 
p < P4, a negative one for p > pa and a maximum at the turning 
point. There exists a limit-incoming geodesic with its turning 
point at the infinite p4 — » 00 for za — » Zi, i-e. with the asymptote 
z-Z\, when p — > 00. 

A test particle with parameters E' , L', and Q' coming from 
infinity (in practice from the accretion disk) parallel to z = to- 
wards the axis of the black hole corresponds to a geodesic which, 
in the limit p — > 00, has an asymptote defined by z = zi - con- 
stant where zi is the impact parameter Therefore it is a limit- 
incoming particle with z < Zi, z < 0, p < 0, ei - 62 - I. In the 
limit p » z, dz/dp - 0, so the tangent has to be parallel to the p 
axis and (IZTT i produces 



Zl 



-I- 

M\E'^ - 1 



1/2 



(52) 



We have plotted in Fig. [3] (see section |5]l an example of a 
geodesic of an incoming particle. We see that, unlike pi , zi does 
not depend upon the black hole parameter a. The incoming parti- 
cles come from the accretion disk, which means that their energy 
E' is rarely very big, which means rarely as big as a or M, which 
characterize the black hole energy, but instead of the same order 
as 1 or VS'. When £' — > 1 then zi — > 00 if (3' 0. However, if 
E' is very big compared to 1 and y/O', then zi — > 0. 

For a given Q', the most energetic incoming particles are 
those with a small impact parameter zi , near to zero. Hence only 
a thin slice of the accretion disk can participate with the greatest 
efficiency in producing Penrose processes leading to the most 
intense possible jet. The point where the ergosphere surface in- 
tersects the z axis is Ze - VA. The value of Ze, for the incoming 
particles, does not play a role like that of p^ for the outgoing 
particles (compare ( l40b to (l52b). 

5. Penrose process and plotting of geodesies 

To make a jet using the geodesies discussed above, we would 
have to assume that incoming particles arrive in the ergosphere 
and undergo a Penrose process. As mentioned earlier, in its origi- 
nal version (Penro se, 69)| each particle may be decomposed into 
two subparticles and one of them may cross the horizon and fall 
irreversibly into the BH, while the other is ejected to the exterior 
of the ergosphere; or the incoming particle may collide with an- 
other particle resulting in one plunging into the BH and the other 
being ejected to the exterior. The second case can correspond to 
a creation of particles, say and from an incoming photon 
(6 - 0) interacting with another inside the ergosphere. We do not 
present here all the possible cases, which are exhaustively stud- 
ied, especially for AGN, in |(Winiams,"95)l|(Wimams, 04)| There 
is also observational evidence for a close correlation between the 
disappearance of the unstable inner accretion disk and some sub- 
sequent ejections from microquasars such as GRS 1915-H1105 
(Mirabel & Rodriguez, 94)1 [(Mirabel & Rodn guez, 99)1 which 
from our point of view could correspond to the instability caus- 
ing disk material to fall through the ergosphere and to then 



give rise to a burst of ejecta from Penrose-like processes. Here 
we are mainly interested in the outgoing particles which fol- 
low geodesies that tend asymptotically towards a parallel to the 
z axis, as described in the earlier section [3] These events are 
closely dependent on the possibilities allowed by the conserva- 
tion equations. In the case when the incoming particle splits into 
two (Ree s et al., 76)| the conservation equations of the energy 
and angular momentum are 



Ein — Equi + Ej^ii, 



(53) 
(54) 



We know from ( l29b that for the outgoing particles we study 
Lz,oiii - 0. The particles falling irreversibly into the BH have 
energy and angular momentum L,jaii. The energy and the 
angular momentum of a particle in the ergoregion can be nega- 
tive which is the basis of the Penrose process. The particle falling 
into the black hole has negative energy, and hence the outgoing 
particle, leaving the ergosphere, has a bigger energy than the in- 
coming particle, 

Eout - Ein + \Efall\- (55) 

We have plotted numerically the geodesies for incoming, 
outgoing and falling particles with the following values for the 
parameters: ajM - 1/2, £■,„ = 200, Eout - 202, Ejaii - -2, 



"Jail 



-100 and Q - Q' = 10 . These values are cho- 



sen in such a way that the geodesies meet inside the ergoregion 
situated in the interior of the ergosphere 



■Pe 



1-1^ 

Pe 



1-1^ 

Pe 



(56) 



and they produce for the asymptotes of the outgoing particles 
pi - 1.64341 and for the incoming zi = 1.58116. The curves 
are built with the initial conditions z[1.55] = 20 for outgoing 
particles and z[0.25] = 0.60 for incoming particles. Their inter- 
section is situated at the point p, = 0.23 and z; - 0.59 inside 
the ergosphere which we can take as the initial condition for the 
falling particle and from where we trace the three curves (see fig. 
4). 




Fig. 3. Penrose process. Plots of the ingoing particle (dashed 
line) coming asymptotically, for p — > 00, from zi = 1.58116, 
to the turning point (p = 0.23, z = 0.59) located inside the er- 
gosphere, whose boundary is indicated, where it decays into an 
infalling particle and an outgoing particle (bold line) following 
a geodesic asymptotic to p = pi - 1.64341 when z — > 00. The 
parameters of the black hole are M - 1 and a - 0.5. The parame- 
ters of the three geodesies are £,•„ - 200, Eout - 202, Efaii - -2, 
Lzjn = Lzjaii = -100, L,,„„, = 0, (3 = (3' = 100000. 
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The exhibition of these numerical solutions with an outgoing 
geodesic which leaves the ergosphere after the Penrose process 
and has vertical asymptote with the value p\ precisely equal to 
( l40b confirms that a model based on such geodesies is possible. 

6. Implications for jet formation 

We have shown that to obtain a jet of particles close to the ro- 
tation axis, it must be formed from particles with (almost) zero 
angular momentum, Ly - 0. If we consider only particles with 
= 0, there is among them a subset which give a perfectly col- 
limated jet, i.e. a set of geodesies exactly parallel to the axis: for 
each allowed value of Q they form a ring of radius 



Pi ^Pe 



1 + 



Q 



aHE^ - 1) 



1/2 



(57) 



whwre - a/M. Note that p = p^ is the circumference of 
the ergosphere in the equatorial plane. For large E the set of all 
these geodesies will give a jet of radius of the order of , with a 
density of particles dependent on the distribution in Q for given 
E, which remains perfectly coUimated all the way to infinity. 

We note that all other geodesies with L- - will spread out 
from the axis along lines z = Kp. An astrophysical jet will of 
course be of only finite extent and not perfectly collimated, so 
it could include such geodesies for suitably large K, as well as 
geodesies with a small L- + 0. 

Thus forming a collimated jet of particles from a Penrose- 
like process, this jet having a narrow opening angle, for a ro- 
tating black hole without an electromagnetic field, depends on 
the initial distribution of particles leaving the ergosphere, or of 
some non-gravitational collimating force, even if we consider 
only particles with L, - 0. 

On the other hand, outgoing particles with small energies, 
namely of the order of their rest energy, £ » 1, and Q > Q have 
asymptotes parallel to the z axis with pi » pg. 

This predicted scale of the region of confined highly ener- 
getic particles might provide a test if the accretion disk param- 
eters provided values for the BH mass and angular momentum, 
in a manner such as discussed in |(McClintock et al., 06 )| and pa- 
pers cited therein, and if the transverse linear scale of the jet near 
the BH could be measured. (Particles of equally high energy may 
exist in p > pi but will spread out away from the axis.) 

Let us make a brief qualitative remark about the observability 
of the two species, (a) and (b), of geodesies outgoing from the 
ergosphere, studied in section 3 (after (45)). As illustrated by 
the figure 2, for each fixed value of pi there is one (b)-geodesic 
only, which is the limit of many (one infinity of) (a)-geodesics 
when the turning point tends to the infinity (z2 °°, Pi ^ Pi)- 
However, the (a)-type geodesies, though much more numerous 
than the (b)-type geodesies, are, directly or indirectly (i.e. by 
radiation, if charged), much more difficult to observe. 

Indeed, contrary to the set of (b)-particles framing the jet 
in one direction (collimation along the poles), the (a)-particles 
ejected from the ergosphere along unbound geodesies at lower 
latitudes are dispersed into the whole 3D-space {An steradians). 
The (a)-particles never produce a beam into one privileged di- 
rection but instead dilute in the whole space. Observed from the 
infinity in one line of sight {9 - constant, (p ^constant), one 
single (a)-particle could directly be detected. While, from the in- 
finity in the line of sight z {9 - Q, V0), the observer will see 
one infinity (each point of the perimeter of the circle of radius 
pi) of (b)-particles. The result is reinforced when we extend it 
to all the possible values of pi . Encircling the foot of the (b)-jet. 



the (a)-particles frame a gerb, from the basis of which a possible 
indirect effect of isotropic radiation emission (from accelerated 
charged particles) could be observed, during the jet eruption. 

Besides, by their dispersion, the pressure the (a)-particles 
locally exert on the ambient medium is much weaker than the 
pressure exerted by the numerous coherent (b)-particles of the 
jet (a narrow parallel beam is more incisive). The (a)-particles 
are probably more rapidly thermalised than the (b)-particles of 
the jet. So, one might expect that many particles ejected at lower 
latitudes never attain infinity (neither the height of the jet), and 
most of them feed the medium, framing a halo around the BH, 
falling inside again, or returning to the accretion disk. 

We noted also that geodesies with <3 < can be asymptotic 
to lines with constant 9. These asymptotes allow us to define 
another type of jet which is bigger and less collimated than the 
previous one. It is interesting to remark that recent observations 
|(Sheth et al., 03)| (Sauty et al., 02) suggest the existence of two 
different types of jets precisely of these sorts, i.e. narrowly and 
broadly collimated. 

There exists an ensemble of geodesies that tend asymptoti- 
cally to these conical characteristics. The unbounded geodesies 
have mainly been discussed, however, by using Boyer-Lindquist 
coordinates r and 9 by the majority of authors. If we rewrite our 
results, using these coordinates, we may interpret our results and 
compare to those of other authors. However, as we show below, 
these coordinates are not as well-suited to the issues we have 
discussed. 

Geodesies with Q may reach low values of pi z, if b2 
is large enough, but must be bounded away from | cos 9\ - I 
(i.e. 9 - or 9 - n), since those values would imply (f- < 0, 
from OJ, cf. KChandrasekhar, 1983)1 p. 348. In practice this 
means that a narrow jet along the axis must be composed of 
particles with very small L.. Particles with non-zero L- could 
only lie within a jet with bounded p for a limited distance, be- 
cause large enough z would imply 9^ < 0. If Q > 0, the or- 
bits reverse the sign of 9 and reach the equatorial plane, and 
would thus be expected to be absorbed by the accretion disk. 
For Q < Q they are confined to a band of values of 9 given by 
the roots of - 0. These are the 'vortical' trajectories of de 
Fehce et al. |(de Felice & Calvani, 72)] |(de Fehce & Curir, 92)| 
|(de Felice & Carlotto, 97) Depending on the maximum open- 
ing angle 9, these may still hit, and presumably be ab- 
sorbed by, a thick accretion disk |(de Felice & Curir, 92)| Such 
orbits can be adequately populated by Penrose-like pro- 
cesses [(Williams, 95)| [(Williams, 04)| and might undergo pro- 
cesses which reduce the opening angle [(de Felice & Curir, 92)[ 
(de Felice & Carlotto, 97)[ A jet composed of such particles 
would tend to be hollow and would have a larger radius p at 
large z than is obtained for orbits with L, - 0, and hence be 
observationally distinguishable. The presence of these escaping 
trajectories spiralling round the polar axis can be associated with 
the gravitomagnetic effects due to the rotation of the hole, one 
of whose consequences is that even curves with L, - have a 
non-zero d<p/dt at finite distances. 

Thus although an infinitely extended jet of bounded p radius 
would only contain particles with L, = 0, which we would ex- 
pect to be a set of measure zero among all particles ejected, we 
shall consider this as a good model even for real jets. In prac- 
tice, interactions with other forces and objects, which would af- 
fect the jet both by gravitational and other forces, have to be 
taken into account once the jet is well away from the BH, and 
these influences might or might not improve the collimation. 
In [(de Fehce & Carlotto, 97)] the authors discussed possible im- 
proved collimation for particles of low L- using forces which 
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have a timescale long compared with the dynamical timescale of 
the geodesies, and which act to move particles to new geodesies 
with changed parameters. It should be noted that if the object 
producing the jet is modelled as a rotating black hole, produc- 
tion of a collimated jet only arises naturally if the object throws 
out energetic particles with low L,, since our discussion shows 
that other particles cannot join such a jet unless there is some 
other strong collimating influence away from the BH. 

However previous authors have not pointed out the existence 
of asymptotes p - pi, presumably because they are less obvious 
when using coordinates r and 9. In fact, considering z — » oo, the 
expressions ( fTSl l and (fTSl l produce cos 9^1- {p^jlz^) + 0{z^^), 
sine ^ (pi/z) + and r a! z + 1 + {p\l2z^) + 0{z-^). With 

these expansions it is clear that in the limit 9 - Q one would have 
to take the limit of r sin 9 to allow p\ to be determined. 

In the same vein, to find the values of asymptotes z - Z\ + 
near the equatorial plane 9 - njl for the incoming particles 
(see (l52T i) one has to study r cos 9 if one uses the coordinates r 
and 9. In fact, one finds for the asymptotic expansion p — > oo 
the following expressions, cosfl ^ (zi/p) + 0(p~^), sin^ =s 1 - 
(z2/2p2) + (9(p-3) and r p+ 1 + + 1 - (a/Mf]/2p + 0(p-\ 



7. Conclusion 

Our main results are the following. 

There are projections of geodesies all over the meridional 
planes. Among these geodesies there are some, with vertical 
asymptotes parallel to z which can form a perfectly collimated 
jet. There are, as well, geodesies with horizontal asymptotes par- 
allel to the radial coordinate p, that can represent the paths of 
incoming particles leaving the accretion disk. 

These two types of geodesies have intersection points that 
can be situated inside the ergosphere. At these points a Penrose 
process can take place, producing the ejection of particles along 
the axis with bigger energies than the energies of incoming parti- 
cles close to the equatorial plane. The energies of outgoing par- 
ticles are significantly larger than the ones of the incident par- 
ticles for the asymptotically vertical geodesies near the scale 
a/M of the ergosphere diameter in the coordinate p, so such 
particles can show collimation around the surface of a tube of 
diameter 2a /M centred on the axis of symmetry. Such colli- 
mated outgoing particles have to have a zero orbital momentum 

= 0, which implies, from the Penrose process, that the in- 
coming particles have a negative orbital momentum, L' < 0. 
Thus the jet has to be fed from incoming particles with ret- 
rograde orbits in the accretion disk. There is evidence for the 
existence of substantial counterrotating parts of accretion disks 
|(Koide et al., 00)| |(Thakar et al., 1997)1 such counterrota- 
tions could explain the viscosity inducing the instabilities which 
trigger the falling of matter towards the ergosphere. It is now 



known |(MirabeI & Rodriguez, 94)| |(MrrabeI & Rodriguez, 99) 
[(Mirabel, 06)| that there is a close connection between instabili- 
ties in the accretion disk and the genesis of jets for quasars and 
microquasars. 

The most energetic incoming particles are those near the 
equatorial plane. Hence the incoming particles which produce 
the most energetic outgoing particles by a Penrose process in the 
ergosphere, whose maximum size is Ze - VA, are those with 
angular momentum L' < and a very small impact parameter 

Also, the limiting diameter of the core of a perfectly colli- 
mated jet depends upon the size of the ergosphere. The effective 



thickness of this part of the jet in this case is of the order of 

2pe = 2fl/M. 

Our idealised model is based on the well-behaved vacuum 
stationary exact solution of Einstein's equations with axial sym- 
metry, namely the Kerr metrics, which does not take into ac- 
count the ambient medium. Though this medium is very dilute, 
it plays a non-negligible role on the more complex global sce- 
nario for jets like progressive widening of the beam, advent of 
knots, lobes, etc... However, for the scenario that we are here 
concerned, namely the beginning of the jet (parsec scale for mi- 
croquasars, while some hundred parsecs for AGN, depending on 
the BH mass), where it is strongly collimated, our approximation 
of test-particles along geodesies is relevant. Indeed, the observed 
jets stemming from active galactic nuclei ejected along the po- 



lar axis have ultrarelativistic speeds, typically 



0.99995c. 



The ejected particles, forming the jets, are thermalized with tem- 
peratures of the order lO^K ( |(Filloux, 09)[ ) producing a lateral 
force from the pressure gradient between the thermal energy 
of the particles in the outflow and the low density enveloping 
medium ((Punsly, 99a) (Punsly, 99b)). The internal particle tra- 
jectories to these jets expand laterally at the speed of sound, be- 
ing of the order = 30km/s [("Filloux, 09)1 asymptotically form- 
ing a conical shape with an opening angle of the order of the 
inverse Mach number (vjvj) - IQ^* radians. As we can see 
((Punsly, 99b), Appendix), the more realistic trajectories corre- 
sponding to such corrective terms represent only a small pertur- 
bation to the geodesies. 

The model that we present to explain the formation and col- 
limation of jets arises essentially from relativistic strong gravi- 
tational field phenomena without resort to electromagnetic phe- 
nomena. From this point of view the model could be interesting 
also for understanding observational evidence of neutral parti- 
cles emitted from the inner jet itself. For example, the recent 
observations of Ultra High Energy Cosmic Rays (diflicult to ex- 
plain, implying neutral particles such as neutrinos, or H or Fe 
atoms, etc. '(Auger, 20 07a)| [(Auger, 2007b)l |(Dermer et al., 08)| 
and (HESS collab., 09)| and references therein) is a new chal- 
lenge. To explain the Very High Energy of such neutral (massive) 
particles, especially neutrinos ,(Auger, 2007a) (Auger, 2007b) 
(which are able to travel freely over large distances), our model 
very naturally suggests that they could be directly coming 
from the collimated inner jet, which would privilege sources 
(BH) with rotational z-axis along the line of sight of the ob- 
servation. Massless particles (photons) would be emitted by 
charged particles accelerated along the collimated inner jet 
KDermer et al, 08)[ [(HESS coUab., 09) ' which would privilege 
sources (BH) with rotational z-axis perpendicular to the line of 
sight of the observation. 

Our model is sufficiently general to fit various types of ob- 
served jets, like GRB, jets ejected from AGN or from micro- 
quasars, whenever they are energetic enough to be explained by 
just a rotating black hole fed by an accretion disk in an axisym- 
metric configuration. The main drawback is the need to prefer- 
entially populate the geodesies which can form such collimated 
jets. Work is in progress on this question to determine a possible 
confrontation of the model with observations. Our preliminary 
studies led us to understand the fundamental role of the function 
P(r) of the geodesies equations (See (2) and (19)). As an exam- 
ple, in the special case where the equation P{r) = has a real 
double root, there exist only two narrow ranges of pi values for 
large values of E. In this case, we can evaluate from the power, 
for example of radio loud extragalactic jets ( [(Willott et al., 99)} , 
or of microquasars jets ( [(Fender et al., 04)) , the particle density, 
the mean kinetic energy by particle, the mean velocity and the 
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Lorentz factor of the jets. These results, since they require a long 
presentation, deserve a separate paper which is under prepara- 
tion. 

The existence of vacuum solutions of the Einstein equa- 
tions of Kerr type but with a richer, not connected, topo- 
logical configuration of the ergosphere (see |(Gariel et al., 02)[ 
figs. 7 to 10), allows us to propose the existence of double 
jets, because they are expected to come out from the ergo- 
sphere. These bipolar jets have been observed (see for instance 
[(Skinner et al797) 1 (Sahai et al., 98)| fig. l, |(Fargion, 03)] fig. 2, 
and |(Kwok et al.,'9 8)) and could be naturally interpreted in a 
generalization of our model. 
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